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ABSTRACT 

Using Chandra X-ray observations of nine nearby, X-ray luminous elliptical galaxies 
with good optical velocity dispersion measurements, we show that a tight correlation 
I exists between the Bondi accretion rates calculated from the observed gas tempera- 

^ . ture and density profiles and estimated black hole masses, and the power emerging 

' from these systems in relativistic jets. The jet powers, which are inferred from the 

, energies and timescales required to inflate cavities observed in the surrounding X-ray 

' emitting gas, can be related to the accretion rates using a power law model of the 

(N ; form log(PBondi/10''^ergs-i ) = A -h Slog (Pjet/lO''^ erg s^^ ), with A = 0.65 ± 0.16 

, and B = 0.77 ± 0.20. Our results show that a significant fraction of the energy asso- 

dated with the rest mass of material entering the Bondi accretion radius {2.2'^q'j per 
cent, for Pjct — lO'^'^ergs"^ ) eventually emerges in the relativistic jets. The data also 
hint that this fraction may rise slightly with increasing jet power. Our results have 
Qh' significant implications for studies of accretion, jet formation and galaxy formation. 

Q . The observed tight correlation suggests that the Bondi formulae provide a reasonable 

■ description of the accretion process in these systems, despite the likely presence of 

^ I magnetic pressure and angular momentum in the accreting gas. The similarity of the 

. Psondi and Pjot values argues that a significant fraction of the matter entering the 

accretion radius flows down to regions close to the black holes, where the jets are 
presumably formed. The tight correlation between Peondi and Pjct also suggests that 
the accretion flows are approximately stable over timescales of a few million years. 
^ ' Our results show that the black hole 'engines' at the hearts of large elliptical galaxies 

and groups can feed back sufflcient energy to stem cooling and star formation, leading 
naturally to the observed exponential cut off at the bright end of the galaxy luminosity 
function. 

Key words: accretion, accretion disks - black hole physics - galaxies: active - galax- 
ies :jets - X-rays: galaxies 
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1 INTRODUCTION 

The black holes in the hearts of galaxies and galaxy clusters 
are commonly observed to be associated with powerful rel- 
ativistic jets. The mechanism by which such jets form and 
the efficiency with which the energy associated with mate- 
rial entering the accretion radius is converted into jet power 
at much smaller radii remains the subject of much debate. 
Such knowledge is important for understanding the nature 
of the accretion process, galajcy formation and the growth 
of supermassive black holes. 

X-ray studies with Chandra and XMM-Newton, build- 



ing on earlier work with ROSAT, have shown that the black 
holes at the centres of galaxies, groups and clusters inter- 
act strongly with their environments, blowing clear cavi- 
ties or 'bubbles' in the surrounding X-ray emitting gas {e.g. 
Bohringer et al. 1993; Fabian et al. 2000, 2003, 2005; Birzan 
et al. 2004; Forman et al. 2005). Simple arguments based 
on the energy required to inflate the bubbles and the sound 
speed of the surrounding gas allow these bubbles to be used 
as calorimeters, measuring the power released by the jets 
{e.g. Churazov et al. 2002). Importantly, the exquisite spa- 
tial resolution of Chandra means that for the largest, nearby 
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ellipticals we are also able to resolve regions close to the 
Bondi accretion radii in these systems. Given estimates for 
the central black hole masses from measured galaxy velocity 
dispersions, the observed density and temperature profiles of 
the X-ray emitting gas can be used to estimate the accre- 
tion rates onto the black holes (Bondi 1952). Taken together, 
the accretion rates and jet powers then allow us to estimate 
the efficiency with which the rest mass of accreted matter is 
converted into jet power. 

The observed luminosities of AGN in elliptical galaxies 
are typically several orders of magnitude lower than would 
be predicted from the standard Bondi formulae, assuming a 
conversion efficiency for matter into radiation of order 10 per 
cent. However, much of the predicted power may actually 
be released in the form of jets. This issue was investigated 
by Di Matteo et al. (2003) and Taylor et al. (2006) who 
noted that although the central AGN in NGC4486 (M87), 
NGC4696 and NGC6166 have bolometric luminosities three 
orders of magnitude below the Bondi predictions (for 10 per 
cent efficiency), the jet powers inferred from the observed 
X-ray cavities are comparable to the Bondi values, within 
a factor of a few. Pellegrini et al. (2003) used related argu- 
ments to infer a minimum conversion efficiency for accreted 
mass into jet power of ~ 1 per cent for the elliptical galaxy 
IC4296. Here, we extend the work of Di Matteo et al. 2003 
and Taylor et al. (2006) to a sample of nine nearby. X-ray 
bright ellipticals with good X-ray data and optical velocity 
dispersion measurements. 

The origin of the observed galaxy luminosity function 
and, in particular, the mechanism responsible for truncat- 
ing star formation in the largest ellipticals is also a sub- 
ject of much debate [e.g. Benson et al. 2003, Croton et al. 
2005; Bower et al. 2005, and references therein). Here, again, 
significant insight has been provided by XMM-Newton and 
Chandra observations {e.g. Tamura et al. 2001; Allen et al. 
2001; Peterson et al. 2001, 2003) which show that the cooling 
rates in the cores of galaxies, groups and clusters are much 
smaller than the values predicted by simple models that do 
not include feedback. Some widespread form of heating must 
be operating in galaxies and clusters to stem cooling. We 
show here that the kinetic power released in jets from the 
central black hole engines can, in principle, provide this heat 
and stem further star formation. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Target selection 

Our target galaxies are bright, nearby ellipticals with accu- 
rate central black hole mass measurements or velocity dis- 
persions which allow the black hole masses to be estimated 
e.g. using the relation of Trcmainc et al. (2002). The Chan- 
dra archive was searched to identify all such systems with 
X-ray luminous gaseous halos and obvious jet-induced cav- 
ities in the X-ray emitting gas. We have restricted our tar- 
get list to those systems in which the bubbles are attached 
(or very close to) to the central AGN, which are therefore 
likely to be undergoing inflation. We have also required that 
the Chandra exposure times be sufficient to allow precise 
measurements of the central density and temperature of the 
X-ray emitting gas. Finally, we required that the X-ray data 



Table 1. Summary of the Chandra observations. Columns list 
the target name (and alternative), distance in Mpc, observation 
date and net exposure after all cleaning and screening procedures 
were carried out. For the Virgo Cluster ellipticals, a distance of 
17 Mpc has been assumed. For NGC 507, 4696, 5846 and 6166 
the distances are luminosity distances for a flat ACDM cosmology 



with Hm = 0.3 and 


a Hubble Constant ifo = 70kms~i Mpc~i. 




cLl (Mpc) 


Date 


Exposure (ks) 


NGC507 


71.4 


2000 Oct 11 


19.1 


NGC4374 (M84) 


17 


2000 May 19 


28.1 


NGC4472 


17 


2001 Jun 12 


21.9 


NGC4486 (M87) 


17 


2000 Jul 29 


33.7 


NGC4552 (M89) 


17 


2001 Apr 22 


53.4 


NGC4636 


17 


2000 Jan 26 


43.2 


NGC4696 


44.9 


2004 Apr 01 


86.7 


NGC5846 


24.6 


2000 May 24 


19.5 


NGC6166 


135.5 


1999 Dec 11 


15.0 



allow us to measure the properties of the X-ray emitting 
gas within one order of magnitude of the Bondi radius, per- 
mitting reliable measurements or extrapolations of the gas 
properties at that radius to be made. These selection crite- 
ria are matched by eight galaxies with data in the Chandra 
archive. For comparison purposes, we also include results for 
one more distant system, NGC6166, the dominant galaxy of 
Abell 2199. Although, due to its greater distance, we do not 
come close to resolving the Bondi radius in this system, the 
density profile follows a simple power-law form which allows 
an interesting extrapolation to smaller radii. 

2.2 Chandra analysis 

The Chandra observations were carried out using the Ad- 
vanced CCD Imaging Spectrometer (ACIS) between 1999 
December and 2004 April. The ACIS-S array was used as 
the primary detector as it offers the best sensitivity to soft 
X-ray emission from the galaxies. The standard level-1 event 
lists produced by the Chandra pipeline processing were re- 
processed using the CIAO (version 3.2.1) software pack- 
age, including the latest gain maps and calibration prod- 
ucts. Bad pixels were removed and standard grade selec- 
tions applied. Where possible, extra information available 
in VFAINT mode was used to improve the rejection of cos- 
mic ray events. The data were cleaned to remove periods of 
anomalously high background using the recommended en- 
ergy ranges and time binning. The net exposure times af- 
ter cleaning are summarized in Table 1. The analyses of 
NGC4486 (M87) and NGC6166 have been reported sepa- 
rately by Di Matteo et al. (2001, 2003). For these objects, 
we have used their density and temperature profiles, scaled 
to the distances summarized in Table 1. 

The temperature and density profiles of the X-ray emit- 
ting gas were determined using the methods described by 
Allen et al. (2004 and references therein). In brief, concen- 
tric annular spectra were extracted from the cleaned event 
lists, centred on the AGN (identified from radio data and/or 
the presence of an obvious X-ray point source in the Chan- 
dra data). These spectra were analysed using XSPEC (ver- 
sion 11.3: Arnaud 1996), the MEKAL plasma emission code 
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Table 2. Summary of the measured Bondi accretion rates and powers. Columns list the galaxy velocity dispersion {a in units kms~^), 
the logarithm of the accretion radius {va in units of parsecs), slope of the inner electron density profile (a), electron density at the 
accretion radius {nc{r = ta) in units cm~^), temperature at the accretion radius {kT{r = ta) in units of keV), logarithm of the Bondi 
accretion rate (Meondi in units of MQyr^^), and the logarithm of the predicted accretion power (PBondi in units of lO^^ergs^^). 
Uncertainties arc 68 per cent confidence limits. Additional systematic uncertainties of 0.2,3 dcx in logrA, and 0.46 dex in logMeondi Eind 
logPeondi) originating from the observed systematic scatter in the logMen— logf relation (Tremaine et al. 2002), may be associated 
with these values. See text for details. 



logo- 



logrA 



a 



nc(r = rA) kT{r = ta) logMeondi logPsondi 



NGC507 
NGC4374 (M84) 

NGC4472 
NGC4486 (M87) 
NGC4552 (M89) 
NGC4636 
NGC4696 
NGC5846 
NGG6166 



2.498 ± 0.014 
2.473 ± 0.010 
2.488 ± 0.013 
2.557 ±0.013 
2.410 ± 0.031 
2.318 ±0.017 
2.418 ±0.014 
2.416 ± 0.012 
2.496 ± 0.012 



1.58: 
1.49 



().09 
0.08 
-0.07 

i.OD_y J2 

^+0.13 
-0.20 

1 26+°" 
Qgg + 0.12 

1 99 + 0.07 
-■-•^^-0.07 

-^•^="-0.08 
,+0.31 
0.20 



2.09^ 



1.31" 



-1.10 
-0.55 
-0.36 

-o.oo; 

-0.83 

-0.31 
-0.63 
-0.35 
-0.50 



0.04 
0.04 
0.19 
-0.18 
+ 0.12 
0.12 
+0.10 
0.10 
+0.08 
0.08 
+0.09 
-0.09 
0.30 
0.30 
0.26 
0.21 
+ 0.07 
0.06 



3.23 
0.94 
0.83 
0.17 
1.30 
0.42 
1.61 
0.35 
0.92 



+ 1.03 
0.71 
0.84 
-0.46 
+0.42 
0.28 
+0.04 
0.03 
+0.76 
0.45 
+0.22 
0.15 
+3.82 
1.15 
+0.52 
0.22 
+0.41 
0.33 



0.74 ± 0.04 
0.71 ± 0.05 
0.70 ±0.16 
0.80 ±0.01 
0.67 ± 0.09 
0.54 ±0.11 
0.81 ±0.05 
0.67 ±0.09 
1.3 ±0.7 



-1.34^ 



1+0.09 
-0.09 

_2 07+0.30 



-1.96 
-1.59 
-2.38 
-3.46 



0.29 
+0.25 
0.23 
+0.28 
0.40 
+0.22 
0.21 
+0.24 
-0.24 
0.56 



-2.35lo.55 

-2 90+''-''^ 
^•»IJ_0.40 
r, 97+0.34 
-0.26 



1.41 
0.69 

0.79 
1.16 
0.37 

-0.71 
0.40 

-0.15 
0.49 



+0.09 
0.09 
+0.30 
0.29 
+0.25 
0.23 
+0.28 
0.40 
+0.22 
0.21 
+0.24 
-0.24 
+0.56 
0.55 
+0.43 
0.40 
+0.34 
0.26 



(Kaastra & Mcwc 1993; incorporating the Fc-L calculations 
of Licdhal, Ostcrhcld & Goldstein 1995) and the photoelec- 
tric absorption models of Balucinska-Church & McCaninion 
(1992). We have included standard correction factors to ac- 
count for time-dependent contamination along the instru- 
ment light path and have incorporated a small correction to 
the High Resolution Mirror Assembly model in CIAO 3.2.1, 
which takes the form of an 'inverse' edge with an energy, 
E=2.08keV and optical depth r = —0.1 (Herman Marshall, 
private communication). Only data in the 0.6 — 8.0 keV en- 
ergy range were used for our analysis. The spectra for all 
annuli for a particular galaxy were modelled simultaneously 
in order to determine the deprojected X-ray gas temperature 
profiles, under the assumption of spherical symmetry. ^ . The 
emission from each spherical shell was modelled as a single 
phase plasma, with the abundances of all metals in each 
shell assumed to vary with a common ratio, Z, with respect 
to Solar values. (The exception to the latter assumption was 
the analysis of NGC4486, reported separately by Di Matteo 
et al. (2003), for which the abundances of key elements were 
permitted to vary independently in each shell.) 

Background spectra were extracted from the blank-field 
data sets available from the Chandra X-ray Center. These 
were cleaned in an identical manner to the target obser- 
vations. In each case, the normalizations of the background 
files were scaled to match the count rates in the target obser- 
vations measured in the 9.5-12keV band. Separate photon- 
weighted response matrices and effective area files were con- 
structed for each region using calibration files appropriate 
for the period of observation. 

Azimuthally-averaged surface brightness profiles were 
constructed from background subtracted, flat-fielded images 
with a 0.492 x 0.492 arcsec^ pixel scale (1x1 raw detector 

^ Note that for NGC4696, we have only analysed data from a 
35 degree wide slice, covering position angles 15-50 degrees. This 
excludes the complex structure to the west of the central source 
discussed by Taylor et al. (2006). For all other sources, the full 
360 degrees of position angle were used. Note that for NGC4486, 
the bright non-thermal emission from the jets was also excluded 
(Di Matteo et al. 2003) 



pixels). Together with the deprojected spectral temperature 
profiles, these were used to determine the X-ray gas density 
profiles (sec e.g. White, Jones & Forman 1997 for details of 
the technique). 



3 CALCULATION OF ACCRETION RATES 
AND JET POWERS 

3.1 Black hole masses and accretion radii 

For the case of M87, the measurement of the central black 
hole mass, Mbh = 3.0 ± 1.0 X 10^ Mq, is from Tremaine et 
al. (2002) and is based on the data of Harms et al. (1994) 
and Macchetto et al. (1997). For the other galaxies in the 
sample, the black hole masses have been estimated using the 
correlation between Mbh and velocity dispersion, cr, given 
by Tremaine et al. (2002) 

log(MBH/M0) =a + /3(logcT/200kms-^ Mpc"^), (1) 

with a = 8.13±0.06 and p = 4.02±0.32. Our analysis uses a 
Monte Carlo algorithm to account for the uncertainties in cr 
and the slope of the log Mbh— log cr relation. Velocity disper- 
sions for the galaxies are drawn from the work of Bernardi 
et al. (2002) and are summarized in Table 2. 

The accretion radius, ta, is the radius within which the 
gravitational potential of the central black hole dominates 
over the thermal energy of the surrounding X-ray emitting 
gas 

TA = 2GMbh/cI (2) 

Here G is the gravitational constant, Cs = \J ^\kT j /imp is 
the adiabatic sound speed of the gas at the accretion radius, 
T is the gas temperature at that radius, /it = 0.62 is the mean 
atomic weight, mp is the proton mass and 71 is the adiabatic 
index of the X-ray emitting gas. The results on the accretion 
radii, summarized in Table 2, account for the uncertainties 
in the central black hole masses (described above) and gas 
temperatures. 

Tremaine et al. (2002) measure an intrinsic (systematic) 
dispersion of 0.23 dex in the log Mbh— log cr relation, which 
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implies a similar systematic uncertainty in the logrA val- 
ues. This systematic uncertainty in dominates over the 
tabulated measurement errors in most cases. The effects of 
accounting for this systematic dispersion are discussed in 
Section 4. 



3.2 X-ray gas properties at the accretion radii 

Fig 1 shows the electron density profiles for the central re- 
gions of the galaxies, determined from the Chandra X-ray 
data. In cases where the X-ray emission from the inner- 
most parts is contaminated by non-thermal emission from 
a central AGN (NGC4374, 4486, 4552, 6166), the affected 
regions have been excluded. The arrows in the figures mark 
the locations of the accretion radii. With the exception of 
NGC 4472, the accretion radii are unresolved by the Chan- 
dra data. However, in general, the inner few points of the 
density profiles can be described by a power-law model of 
the form ne(r) oc r" with -1.1 < a < (Table 2). We have 
used the power-law models to estimate the gas densities at 
the accretion radii, in each case accounting for uncertainties 
in the slopes and normalizations. For M87 (NGC4486), the 
density profile fiattens within r ~ lOrA and so we have as- 
sumed that the density at rx is equal to the value at the 
innermost measurement radius. (We assign a systematic un- 
certainty in the slope a = 0.0 ±0.1 in this case.) Note that 
the regions fitted with the power law density models are 
the regions for which we expect the assumption of spherical 
symmetry to be most valid, based on the Chandra images. 

The temperature profiles for the inner regions of the 
gala:x;ies are shown in Fig. 2. In general, the profiles do not 
exhibit steep gradients within r ~ 10— lOOrA. We have there- 
fore assumed that the temperatures at the accretion radii 
can be estimated from the values at the innermost measure- 
ment radii. The temperature results axe summarized in Tar 
ble 2. Systematic uncertainties in the temperature estimates 
are unlikely to affect other results reported here significantly. 



3.3 Calculation of the accretion rates 

The calculation of accretion rates onto the- central black 
holes follows the work of Bondi (1952). Similar calculations 
are described by e.g. Di Matteo et al. (2001, 2003), Chura- 
zov et al. (2002), Pellegrini et al. (2003) and Taylor et al. 
(2006). 

Under the assumption of spherical symmetry and negli- 
gible angular momentum, the rate of accretion of the X-ray 
emitting gas at the accretion radius can be written as (Bondi 
1952) 

Meondi = 47rA(GMBH)^cJ^p = TrAcspri, (3) 

where p is the density of the gas at the accretion radius 
(we assume p = l.lSnemp) and A is a numerical coefficient 
that depends upon the adiabatic index of the accreting gas. 
For an assumed adiabatic index 71 = 5/3, A = 0.25 (Bondi 
1952). Note that for 71 = 5/3, p/cf is constant within the 
accretion flow, making it appropriate to use the gas temper- 
ature and density measured at the accretion radius in our 
calculations, in place of the values at infinity. (In essence, we 
are modelling the accretion flows as Bondi flows inward of 
the accretion radii. For r > ta, the gravitational potentials 



of the galaxies start to dominate and the Bondi formulae 
become inappropriate.) 

For an efficiency r], relating the accretion rate at the 

accretion radius to the total energy emitted from within 
that radius, the ma^ximum power released from the black 
hole system is 



-pBondi = JjMBondiC 



(4) 



The values of Peondi, for an assumed value of 77 = 0.1 are 
summarized in Table 2. The statistical uncertainties in the 
results have been estimated from Monte Carlo simulations, 
which account for the uncertainties in the a measurements, 
the slope of the logMeH— logcr relation and uncertainties in 
the temperature and density measurements at the accretion 
radii. 

An intrinsic (systematic) dispersion of 0.23 dex in the 
log A'/bh— log CT relation (Tremaine et al. 2002) implies a 
systematic uncertainty of 0.46 dex in both logMsondi and 
logPsondi- These systematic errors are comparable to the 
measurement errors for most galaxies, and dominate for 
NGC507. 



3.4 Calculation of jet power 

To estimate the kinetic power of the jets, we first estimated 
the energy, E required to create the observed bubbles in the 
X-ray emitting gas. For slow expansion rates, this is the sum 
of the internal energy within the bubble and the PAV work 

done 



72 - 1 



PV + PV 



72 



72 ■ 



1 



PV, 



(5) 



where P is the thermal pressure of the surrounding, X-ray 
emitting gas (which can be determined from the observed 
X-ray temperatures and densities), V is the volume of the 
cavity and 72 is the mean adiabatic index of the fiuid within 
the cavity. For the case of a bubble filled with relativistic 
plasma, 72 = 4/3 and E = 4PV. Some additional energy 
will also be transfered into sound waves (Churazov et al. 
2002; Fabian et al. 2003, 2005, 2006; Forman et al. 2003; 
Ruszkowski et al. 2004) . We have neglected this extra energy 
in our calculations, although the power involved is likely to 
be smaller than the quoted uncertainties on the jet powers. 

The regions identified with the bubbles are shown in 
Fig. 3. In general, the bubbles appear to be 'attached' to 
the central AGN and are likely to be undergoing expansion 
at present. A possible exception to this is NGC4472, where 
the bubbles may have recently detached. In cases where the 
edges of the bubbles are not clearly defined in the X-ray 
data, we have used radio observations to estimate their vol- 
umes (Fig. 3). The X-ray and radio data show that in gen- 
eral, the bubbles or cavities are approximately elliptical. We 
have therefore parameterized the projected areas of the bub- 
bles in terms of a semi-axis length, ri, along the jet direction 
and semi-axis width, rw, across it. This allows the cavities 
to be modelled as ellipsoids with volumes V = 47rrirwr'd/3, 
where rd is the unknown depth of the cavity along the line 
of sight. We model the cavities as prolate ellipsoids with 
V = Awnr^/S, although systematic uncertainties in n, rw 
and, especially, rd mean that our estimates of the bubble 
volumes should be regarded as uncertain at the factor ~ 2 
level. 
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Figure 1. The observed electron density profiles for the central regions of the galaxies, as determined from the Chandra data. For NGC 
4374, 4486, 4552 and 6166 the data from the innermost radii, which are contaminated by non-thermal emission from the central AGN, 
have been excluded. The dashed lines show the best- fitting power law models (ne(r) = nir") which have been fitted to the inner few 
bins; the number of bins fitted are noted in parentheses, after the galaxies names below. The Bondi accretion radius, r^., for each galaxy 
is marked with an arrow. From top to bottom, left to right, we show NGC507(10), 4374(3), 4472(4), 4486(1), 4552(3), 4636(4), 4696(4), 
5846(3) and 6166(14). 



Following Birzan et al. (2004; see also Dunn & Fabian 
2004, Dunn, Fabian & Taylor 2005) we have estimated the 
ages of the bubbles as 

tagc = R/Cs, (6) 

where R is the distance of the bubble centre from the 
black hole, which provides a reasonable match to some nu- 
merical simulations. We note that the bubbles have lower 



mass density than the surrounding X-ray emitting gas and 
can therefo re be expe cted to rise with a buoyancy veloc- 
ity, «b = ^j2gVjSCh). Here S = Trr^ is the cross-sectional 
area of the bubble (in the rise direction) , V is the volume, 
Cd ~ 0.75 is the drag coefficient and g = GM{< R))/R^ is 
the gravitational acceleration. The times required for bub- 
bles of the observed sizes to rise buoyantly through a uniform 
medium from the galaxy centres to their current radii is ap- 
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Figure 2. The X-ray gas temperature profiles for the central regions of the galaxies determined from the Chandra data. 



proximately tbuoy = R/vh- Application of this formula leads 
to comparable timescales to the tago values. However, since 
the bubbles studied here are typically still attached to the 
jets, we adopt tagc as our relevant timescale. 

The power involved in 'blowing' the bubbles can be es- 
timated as 

Pjct = iSAagc (7) 

The results on Pjot are summarized in Table 3. The un- 
certainties on Pjct are determined from Monte Carlo simu- 
lations which include the sources of statistical uncertainty 



mentioned above as well as a factor 2 systematic uncertainty 
in the bubble volumes. 



4 THE RELATION BETWEEN ACCRETION 
RATE AND JET POWER 

Fig 4 shows the results on the accretion power Paondi, for 
an assumed efficiency rj = 0.1 for the conversion of accreted 
rest mass into energy, versus the jet power, Pjet, determined 
from the bubble properties. Both powers are in units of 
10'*'^ergs~^ . For each galaxy, we have combined the power 
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Figure 3. Chandra X-ray images for the central regions of the galaxies. The X-ray images are constructed in the 0.5-8 keV band and 
have been adaptively smoothed using 2 or 3-pixel Gaussian kernels. The regions identified with the X-ray cavities, which are used to 
estimate the jet powers, are indicated by bold magenta contours. The narrow black contours show the 1.5 GHz radio emission from 
archival Very Large Array (VLA) data. The radio data help to identify the cavities in cases where identifications from the X-ray data 
alone may be ambiguous. The radio contours are logarithmically spaced and have been adjusted to highlight the extended lobe emission. 
The horizontal scale bar in the bottom left corner of each image denotes a length of Ikpc, except for NGC507 where the scale is lOkpc. 



estimates for both bubbles, with the uncertainties added in 
quadrature. A clear correlation between the Psondi and Pjct 
exists, which can be described by a power law model of the 
form 



log 



(10« ergs"^ 



= A + Bloe 



rjot 



(10« ergs-1 



(8) 



Using the fit statistic, which accounts only for errors 
in the Psondi values, we find A = 0.67 ± 0.07 and B = 
0.74 ± 0.08, with = 5.3 for 7 de grees of freedom (the 
quoted uncertainties on A and B are 68 per cent confidence 
intervals). Using the BCES(F|X) estimator of Akritas & 
Bershady (1996), which accounts for errors in both axes and 
the presence of possible intrinsic scatter, we obtain A = 
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Table 3. Summary of the measured jet powers. In order, we list the distances R from the black holes to the centres of the X-ray bubbles, 
the energy E required to expand the bubbles, the bubble ages tage, and the power required to inflate the bubbles over this timescale 
Pjet • The quoted uncertainties in R include statistical measurement errors only. A factor 2 systematic uncertainty in the bubble volumes 
has been included in the calculation of E and Pjet . 



Galaxy 


Lobe 


R 


E = ApV 


^agc 


-Pjot 






(kpc) 


(10''4 erg) 


(10« yr) 


(10*^ ergs-l ) 




E 


12.8 ± 1.0 


4o9UU ± ziiuU 


OO 1 O T O 

zo.D ± 2.16 


D.60 ± 2.91 




W 


11.3 ± 1.0 


293U0 ± 13300 


25.6 ± 2.18 


3.62 ± 1.70 


T^i /t o ^ A 


N 


o r\n 1 n o 

6.\jA ± U.o 


IDOO ± 693 


6.07 ± 0.56 


n o/^o 1 r\ OT1 

0.868 ± 0.371 




S 


A or' 1 C\ A 
4.0D ± U.4 


2050 ± 854 


10.0 ± 1.08 


0.658 ± 0.268 


JNCjC447z 


E 


/loo 1 r\ A 

4.33 ± 0.4 


Ton 1 OOT 

730 ± 287 


5.28 ± 0.56 


0.445 ± 0.181 




W 


3.40 ± 0.3 


606 ± 244 


5.29 ± 0.57 


0.362 ±0.150 


M87 


w 


1.86 ±0.05 


1210 ± 782 


2.92 ±0.13 


1.26 ±0.84 




E 


1.51 ±0.05 


1670 ± 1100 


2.38 ±0.11 


2.18 ± 1.44 


NGC4552 


N 


0.58 ±0.02 


35.0 ± 12.7 


1.53 ±0.08 


0.0727 ± 0.0274 




S 


0.58 ±0.02 


39.4 ± 14.3 


1.51 ±0.08 


0.0829 ± 0.0307 


NGC4636 


NE 


0.69 ±0.07 


8.10 ± 3.21 


1.77 ±0.15 


0.0147 ±0.0060 




SW 


0.66 ±0.07 


8.04 ±3.25 


1.68 ± 0.17 


0.0152 ±0.0060 


NGC4696 


E 


2.49 ± 0.3 


601 ± 338 


6.06 ± 0.66 


0.313 ±0.177 




W 


3.00 ±0.3 


826 ± 438 


5.54 ± 0.60 


0.478 ± 0.244 


NGC5846 


N 


0.72 ±0.12 


21.0 ± 11.2 


1.81 ± 0.30 


0.0367 ± 0.0185 




S 


0.72 ±0.12 


20.8 ± 11.3 


1.81 ±0.30 


0.0374 ± 0.0188 


NGC6166 


E 


2.30 ±0.10 


829 ± 336 


3.51 ±0.27 


0.741 ±0.316 




W 


2.40 ±0.10 


981 ± 383 


3.66 ± 0.28 


0.838 ± 0.353 



0.64 ± 0.08 and B = 0.78 ± 0.15. The mean deviation about 
the best fitting BCES model is <T(logPBondi) = 0.11. 

The results on A and B quoted above do not ac- 
count for the effects of intrinsic (systematic) scatter in 
the log Mbh— logo- relation. We have used further Monte 
Carlo simulations to examine these effects, introducing an 
intrinsic dispersion of 0.46 dex in logPsondi. Fits to the 
simulated data sets using the BCES(F|X) estimator give 
A = 0.65 ± 0.16 and B = 0.77 ± 0.20, in good agreement 
with the results presented above but with (slightly) larger 
error bars on the fit parameters. 

The results shown in Fig 4 indicate the presence of a 

strong correlation between Psondi and Pjct. A power- law 
model provides a good description of the data. However, 
in gauging the origin and significance of this correlation, 
we must also consider effects that could arise from the 
plotted quantities having factors in common. Both Peondi 
and Pjet depend on the distances to the objects. However, 
Peondi oc (^l"'^, whcrcas Pjct tx d]^^. Thus, the distance de- 
pendences cannot lead to the observed positive correlation. 
Both axes in Fig 4 also depend on the temperature T of 
the X-ray emitting gas. However, the temperature shows 
little variation from object to object and varies only mildly 
with radius between the accretion and bubble radii. More- 
over, the dependences on temperature are approximately 
Peondi oc T~i-5 and Pjet oc T^'^, which cannot lead to the 
observed positive correlation. Finally, both axes involve the 
gas density, measured at the accretion radius for Peondi and 
at the bubble centres for Pjet. These radii are very different 
and the density profiles vary significantly from object to ob- 
ject, meaning that the densities are essentially uncorrelated. 
(In detail, a mild anti-correlation is observed). 



We conclude that the observed correlation between 
ffeondi and Pjet indicates a tight physical connection between 
the two quantities. 



5 IMPLICATIONS AND DISCUSSION 

We have shown that for supermassive black holes at the 
centres of large. X-ray luminous elliptical galaxies, a remark- 
able, tight correlation exists between the 'Bondi' accretion 

rates inferred from the Chandra X-ray data and observed 
galaxy velocity dispersions, and the power emerging from 
these systems in relativistic jets. Our result has important 
implications for the nature of the accretion process and for 
issues relating to feedback, the growth of black holes and 
galaxy formation. 

The relationship between the Bondi accretion power 
and jet power can be described by a power law model of the 
form logPsondi = 0.65(±0.16) + 0.77(±0.20)log Pjct, where 
PBondi = O.lMBondiC"^ and Pjot is the power associated with 
inflating the cavities and providing the internal energy of 
the plasma that fills them. A significant fraction (2.21q;° 
per cent, for Pjet = 10'*^ergs~^) of the energy associated 
with the rest mass of material entering the Bondi accretion 
radius emerges from the systems in relativistic jets. There is 
a slight indication that this fraction increases as Pjet rises, 
from l-Slg g per cent at Pjet = lO**^ ergs~^ , to 3.7^.1 ^ per 
cent at Pjet = lO''* ergs"'^ . (Here, the quoted uncertain- 
ties include all sources of statistical and systematic error 
discussed in the text, including systematic scatter in the 
log Mbh— log CT relation.) 

The existence of such a tight correlation suggests that 
the Bondi formalism provides a reasonable description of 
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Figure 4. The logarithm of the Bondi accretion power (in units 
of 10^^ erg s^-"- ) determined from the Chandra X-ray data, as- 
suming an efhciency »j = 0.1 for the conversion of rest mass into 
energy, versus the logarithm of the jet power (also in units of 
10*^ ergs~^ ). The dashed line shows the best-fitting linear-plus- 
constant model, determined using the BCES estimator. (See text 
for details.) 



the accretion process in these systems, despite the fact that 
the accreting gas contains magnetic fields {e.g. Taylor et al. 
2006) and, presumably, has some angular momentum. More- 
over, the similarity of the Pjct and Psondi values suggests 
that a significant fraction of the matter passing through 
the accretion radius flows all the way down to regions close 
to the black hole, where it presumably provides the power 
source for the jets. This provides an important constraint on 
both accretion and jet formation models. In particular, our 
results limit the amount of material that may be lost from 
the accretion flows en route to the region of jet formation 
and requires that jet formation must be efficient; a few per 
cent of the energy associated with the rest mass of material 
entering the accretion radius emerges in the jets. 

The origin of the jets may be related to the spin of 
the central black hole [e.g. Rees 1978; Hughes & Blandford 
2003). In that case, the observed tight correlation between 
jet power and accretion rate may suggest a narrow range of 
black hole spins for the objects in our sample. The X-ray 
emitting gaseous halos in the giant elliptical galaxies stud- 
ied here are likely to have existed for billions of years. The 
growth of the black holes may therefore have been domi- 
nated by gas accretion, which tends to lead to large spin pa- 
rameters {e.g. Volonteri et al. 2005) and it is possible that 
the black holes are rotating at close to their maximal rates. 
However, the detailed interaction between accretion, spin-up 
and jet power remains to be explored. 

The tight correlation between Psondi and Pjet also sug- 
gests that the accretion flows in the galaxies are stable over 
the periods required for gas to flow from the accretion radii 
to the base of the jets plus the times required for the bubbles 



to be inflated. The flow times tsow ~ ta/cs are of the order 
a few 10'' to a few 10® yr, whereas the bubble ages range 
from 10® - 2 X lO^'yr (Table 3). Thus, our results suggest 
that the accretion flows in the present sample of objects are 
approximately stable over timescales of 10® — lO'^ yr. 

The black holes studied here are accreting at rates of 
a few 10~* to a few lO"'^ yr~^ and converting a few 
per cent of their accreted rest mass into jet energy. The nu- 
clear luminosities (Pellegrini 2005) are typically 2-4 orders 
of magnitude lower than the Pjet values, so the systems are 
'quiescent' in terms of their nuclear luminosities. (Orienta- 
tion and obscuration effects will play some role although are 
unlikely to modify this conclusion e.g. Di Matteo et al. 2001, 
2003; Allen, Di Matteo & Fabian 2000.) It is interesting to 
note that the measured jet powers are comparable to the 
bolometric luminosities predicted by standard, radiatively 
efficient accretion disk models with r; ~ 0.1 (Shakura & 
Sunyaev 1973). Thus, it is possible that the ratio of accre- 
tion rate/total emitted power (bolometric luminosity plus 
jet power) may remain relatively stable as the accretion 
mode varies from a quasar to optically quiescent phase. 

Our results have significant implications for models of 
galaxy formation and support the idea that feedback from 
central black holes is important in shaping the bright end of 
the galaxy luminosity function (see e.g. Croton et al. 2005 
and Bower et al. 2005 for recent discussions). For our sam- 
ple, the observed efficiency with which energy associated 
with the rest mass of accreted matter is fed back into the 
surrounding X-ray emitting gas via the jets is similar to that 
assumed in the semi-analytic models of e.g. Croton et al. 
(2005). Those models have been shown to be able to explain 
the exponential cut off at the bright end of the galaxy lumi- 
nosity function and the fact that the most massive galaxies 
tend to be bulge dominated systems in clusters, containing 
systematically older stars. 

Our results show that the supermassive black holes in 
elliptical galaxies are prodigious energy sources, with the 
bulk of the energy emerging in the form of relativistic jets 
rather than radiation. In order for this energy to be cou- 
pled back into the surrounding matter, the jets must have a 
working surface which is provided in these cases by the X- 
ray emitting gas. For the elliptical galaxies studied here, the 
power fed back into the X-ray emitting gas is, in principle, 
sufficient to offset cooling and heat the gas, thereby prevent- 
ing further star formation. (Note that Best et al. 2006 apply 
simple scaling arguments to a large elliptical galaxy sam- 
ple to show that the time-averaged power output from radio 
sources can, in principle, balance radiative cooling.) Since 
we expect all large galaxies to contain a supermassive black 
hole, even at high redshifts, this suggests that a key stage 
in the formation of such galaxies will be when they first be- 
come large enough to maintain an X-ray gaseous halo. (This 
then provides a working surface for the jets, allowing feed- 
back and preventing further galaxy growth. Although the 
jets could have been providing power before this point, they 
are unlikely to have had an effective working surface.) Note, 
however, that although star formation is truncated when jet 
feedback becomes efficient, the black hole may continue to 
grow via accretion from the X-ray gas at approximately the 
Bondi rate. 

A closely related topic is that of cooling flows in galaxies 
and clusters, and the mechanism responsible for preventing 
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gas cooling to low temperatures {kT < IkeV). As stated 
above, the Pjet values in the galaxies and groups studied 
here are sufficient to offset further cooling of the X-ray gas. 
It is interesting to consider extending the correlation be- 
tween Peondi and Pjet to larger systems where the cooling of 
the gas is more rapid. For the largest cooling flow clusters, 
the cooling luminosity exceeds 10''^ ergs~^ . If the black hole 
masses and accretion rates in these galaxies are not signifi- 
cantly larger than in NGC507 and M87, then the jet power is 
unlikely to be sufficient to balance cooling and some cooling 
to very low temperatures and associated star formation can 
bo expected to occur (sec also Fabian et al. 2002). This is 
consistent with observations indicating the presence of vig- 
orous star formation and large molecular gas concentrations 
in the largest cooling flow clusters, but relatively little in 
less X-ray luminous systems {e.g. McNamara & O'Connell 
1992; Allen 1995; Crawford et al. 1999; Edge 2001). 

Finally, we note that although systematic uncertainties 
have been included in the extrapolation of the gas properties 
to the Bondi radii and in our estimation of black hole masses 
and bubble volumes, other sources of systematic uncertainty 
remain. For example, the assumption that A = 0.25 in the 
calculation of Peondi may not be precise, and our assumption 
that the bubble ages can be estimated as tage is probably 
uncertain at the factor ~ 2 level. However, we do not expect 
such uncertainties to affect the main conclusions drawn here. 
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